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A technique to deduce the virtual temperature from the combined use of the 
equations of fluid dynamics, observed wind and observed radiances is described. 
The wind information could come e.g., from ground-based sensitive very high 
frequency (VHF) Doppler radars and/or from space-borne Doppler lidars. The 
radiometers are also assumed to be either space-borne and/or ground-based. 

From traditional radiometric techniques the vertical structure of the 
temperature can be estimated only crudely. While it has been known for quite 
some time (GAL-CHEN, 1978; HANE and SCOTT, 1978) that the virtual temperature 
could be deduced from wind information only, such techniques had to assume the 
infallibility of certain diagnostic relations. The proposed technique is an 
extension of the Gal-Chen technique. It is assumed that due to modeling 
uncertainties the equations of fluid dynamics are satisfied only in the least 
square sense. The retrieved temperature, however, is constrained to reproduce 
the observed radiances. It is shown that the combined use of the three sources 
of information (wind, radiances and fluid dynamical equations) can result in a 
unique determination of the vertical temperature structure with spatial and 
temporal resolution comparable to that of the observed wind. 


1. INTRODUCTION 

A few years ago, GAL-CHEN (1978) and HANE and SCOTT (197 8) noted that if 
sufficiently accurate measurements of the wind and its time history could be 
obtained from Doppler radars, this information would, in principle, define also 
the thermodynamic structure. In essence, this is done by requiring that the 
data will also satisfy the momentum equations in the least square sense. This 
has reduced the problem to a classical calculus of variation problem (COURANT 
and HILBERT, 1953). The form of the momentum equations assumed in these 
studies is quite general and is, in principle, applicable to small meso- and 
large-scale atmospheric motions. While not immediately obvious, when the 
approximations appropriate to large-scale atmospheric flows are employed, the 
above variational formulation is reduced to solving a classical balance 
equation (HALTINER and WILLIAMS, 1980) of obtaining the geopotential from the 
wind. 


GAL-CHEN and KROPFLI (1984), ROUX et al. (1984) and HANE and RAY (1985) 
have tested the practical utility of the above-mentioned variational formu- 
lation on a variety of observed small-scale phenomena; planetary boundary layer 
(PBL) convection in the Gal-Chen and Kropfli case; severe storms for the Hane 
and Ray case, and a tropical squall line for the Roux et al. case. In all, the 
three case studies of temperature and pressure are deduced from observed 
Doppler radars wind. Satisfactory agreement with in situ thermodynamic 
observations is reported in all three cases. 

As GAGE and BALSLEY (197 8) point out, sensitive Doppler radars can be used 
to obtain mesoscale wind profiles under all weather conditions. The vertical 
resolution is up to 100 m. The time resolution is about 1 hour and the 
horizontal resolution is determined by the average distance between the 
profilers. Comparable resolution is not obtainable from radiometric 
measurements of the atmosphere either from the ground or from satellites. The 
purpose of this paper is to extend and modify the GAL-CHEN (197 8) technique to 
satisfy the following requirements; 
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(a) The horizontal momentum equations are satisfied in the least 
square sense (to be defined further below). 

(b) * The hydioslatic dqhsWaiirfc is satisfied exactly. 

(c) The 4 thfemoaynamic^eq&atiori is satisfied in the least square 
sense. 

(d) The radiative transfer equation at various frequencies is 
satisfied exactly. 

(e) Given wind and radiances as input, virtual temperatures 
should be obtained as an output. The retrieved temperature 
should have horizontal and vertical resolution compared to 
that of the observed wind. 

In this paper only the theory is developed. The practical utility 
remains to be checked. This should be done first by simulation studies and 
then by examining real data. The task is vast and difficult and I hope that 
the theory developed here will stimulate other researchers to check its 
practical utility and to seek even better ways to estimate the virtual 
temperature. 

The technique discussed in this paper has some similarities to the 
techniques considered by KUO and ANTHES (1985) and by BRUMMER et al . (1984). 
However, it also has some potentially important differences. These include 
inter alia : 

(a) Lateral boundary conditions for the temperature are obtained 
directly from the wind data rather than prescribed from a 12- 
hour forecast as in Kuo and Anthes or, as in the B rummer et al. 
case, prescribed from a vertically smoothed temperature profile 
obtained by pure radiometric techniques. 

(b) In both the Kuo and Anthes and Brummer et al . techniques, the 
horizontal divergence equation is used; as a weak constraint 
in the Brummer case and as a diagnostic equation for the geo- 
potential in the Kuo and Anthes case. In our case, attempt 
is made to satisfy, albeit in the least square sense, all the 
prognostic equations relevant to describing mesoscale motions, 

(c) The Kuo and Anthes approach does not utilize the information 
contained in the radiances. It is as sumed that in nature the 
divergence equation is satisfied exactly. This is not true even 
if the wind measurements are error free. In the Brummer 
technique, a temperature profile is sought that will, on the 
one hand, satisfy the divergence equation as close as possible, 
and on the other hand, is also not too far from the smooth 
temperature profile retrieved from radiometric data. Our 
technique, however, demands that the retrieved temperature satisfy 
the radiative transfer equation, augmented by additional dynamical 
constraints. Unlike the pure radiometric techniques, the above 
set is mathematically well posed and no a priori smoothing or 
statistical constraints need to be imposed on the retrieved 
temperature. 

2. MODELING ASSUMPTIONS 

Governing equations . The governing hydrostatic primitive equations in 
Cartesian x, y, z coordinates may be written as: 

Continuity equation, 

Dp/Dt + PV-u = 0 


(i) 
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Horizontal momentum equations, 

Du/Dt = -(l/p)3p/3x + + fv (2) 

Dv/Dt = -(l/p)3p/3y + F 2 - fu (3) 

Here D/Dt is a symbol for total derivative 

D/Dt = 3/3t + u3/3x + v3/3y + w3/3z (4) 

f is the Coriolis parameter (f = 2ftsin<|>; ft is the earth angular velocity and 
4> is the latitude). For convenience we are displaying the equations using 
Cartesian coordinates. Nevertheless, the extension of our ideas to spherical 
coordinates is obvious and all our subsequent discussions (conclusions) are 
valid for spherical geometry. The hydrostatic equation is given vis., 

3p/3z = - pg (5) 

An approximate form of the thermodynamic equation neglecting the contribution 
of moisture to the density p and to the heat capacity (under constant pressure) 

C is 
P 

C p DT/Dt - d/p)Dp/Dt = S h (6) 

The heat capacity under constant pressure is given vis. C = (7/2) R with R 
the gas constant for dry air. The equation for conservation of water vapor is 

Dq/Dt = S v (7) 

The equation of state is 

p = pRT (8) 

r v 

Here u is the three-dimensional wind vector £ = (u, v, w) . u is the horizontal 
velocity in the x direction; v is the velocity in the y direction and w is the 
velocity in the vertical direction, z. The density of dry air is denoted by 
P; p is the pressure; g is the acceleration of gravity; F^ and F« are 
symbols for turbulent friction forces (of dimension Newton/kg) which in this 
study we asstime that they can be either measured directly or parameterized 
based on wind observations. T is a symbol for temperature; p is the pressure. 
T^ is the virtual temperature defined vis. 

T - T = 0.61 qT 

V 

Here q is the water vapor mixing ratio (expressed in 10”3g/(kg of dry air)). 

The symbol is for sources (or sinks) of heat energy. Since in this study 
we are limiting ourselves to relatively short time scales (0-12 hours), 
radiative processes are presumed to be of secondary importance (SMAGORINSRY, 
1974) and the major source of heating in the free atmosphere is due to 
precipitation. The major heating source in the planetary boundary layer (PBL) 
is assumed to be fluxes of sensible heat. is a symbol for sources or 

sinks of water vapor and in accordance with our previous presumption that the 
major contributor to in the free atmosphere is the removal of vapor by 
precipitation. In the PBL, the major source is evaporation from the ground. 

The nature of the data and/or the parameterizations. Our ma j or 
assumptions about the nature of the observed data or the parameterizations 
employed are as follows: 
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(a) Horizontal motions can be measured by means of powerful Doppler radars 
(frequency range is 50-900 MHz), The measurements have accuracy of 

+1 ms~l; are such that all motions with time scales with less 
than 1 hour have been filtered and are possible under all weather 
conditions (LITTLE, 1982). 

(b) Vertical motions with scales described in (a) can be either deduced 
from the horizontal motions (using the mass continuity equation) or 
else can be measured directly by Doppler radars (NASTROM et al., 

1985). To be useful for predictions of synoptic scale motions, the 
accuracy of the deduced (or measured) vertical motions must be of the 
order of + 1 cms~l (HALTINER and WILLIAMS. 1980). 

(c) Remote sensing of temperature and moisture profiles using ground- 
based and/or space platforms renders some useful information under 
almost all weather conditions (WESTWATER et al . , 1985). (This is true 
only if the infrared channels are augmented by additional channels 
from the microwave. Otherwise, contamination from clouds may be 
severe. Furthermore, microwave measurements are contaminated under 
the presence of heavy rain.) The temperature and moisture retrieved 
from these measurements typically have poor vertical resolution. As 
a result, the retrieved profiles have an accuracy of no better than 

+ 3°C for temperature and + 5 g/kg for the moisture. 

(d) At the minimum, it is assumed that the measurements described in (a)- 
(c) are available in at least three spatial locations to be able to 
define a triangle. The satisfaction of this requirement would enable 
calculations of horizontal gradients. It must be borne in mind that 
the distance between the stations also determine the smallest scales 
that can be resolved by such a network. Thus, even though the 
horizontal wind measurements described e.g., in (a) may contain 
spatial scales of motions smaller than the distance between the 
stations, the computed horizontal gradients cannot properly resolve 
this information, 

(e) As is customary in numerical weather prediction (NWP) models 

(HALTINER and WILLIAMS, 1980), we assume that all motions and 
processes with spatial and temporal scales that cannot be resolved 
by the network can either be "parameterized" in terms of what is 
observed or measured directly. For instance, F.^ and F 2 (in (2) 
and (3)) which are turbulent friction terms may be estimated from 
single Doppler radar data (KROPFLI. 1984). Alternatively, one may 
attempt to parameterize it in terms of the larger scale winds. (The 
simplest parameterization is to set ? = f 2 =0.) Another example 

is the precipitation rate and the vertical distribution of latent heat 
release which may be evaluated using conventional radars (DOVIAK, 

1981) or from satellite data (ATLAS and THIELE. 1981) or parameterized 
(e.g., ignored). 

3. ALGORITHM DEVELOPMENTS 


Deduction of horizontal virtual temperature gradients . Taking into 
account our assumption (a) and (b) in the previous section, we may write the 
horizontal momentum equations (2) and (3) as: 


(1/P)V H P = Jg 


Here, G = (Gh,G 2 ) is a given two-dimensional vector 
-Du/Dt + fv + F ; G = -Dv/Dt - fu + F ) which, in p 


(G 1 " 


function 

v» 2 — 4 sv/.ls«- a.- . u 2 / wnjL^ti, in principle, can be 

computed from the observed wind; V is the two-dimensional gradient 
operator. Differentiating (10) with respect to 
and the equation of state (8) one gets 


(l/P) v „P3lnT /3z + gV InT = 3 g/3z. 
n v ri v 


( 10 ) 


z and using hydrostatic (5), 
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Taking it into account (10) 


G31nT v /3z + g^lnTy = (11) 

(11) can be considered as a generalization of the thermal wind relation in z 
coordinates. In fact, for g = (fv;-fu) the thermal wind is reproduced. 


Equation (11) expresses horizontal and vertical temperature gradients in 
terms of observed quantities (i.e., winds and its derivatives). Together with 
other relations to be used further below, it can be used to infer the vertical 
structure. Nevertheless, it is also useful to consider several approximations 
of (11). First consider the ratio (denoted by R ) of gV InT to 
G31nT /3z. Traditional scale analysis considers 3 (e. g. , PEDlXsKY, 1979, pp. 
i'-10) V dictates that the order of magnitude of the above-mentioned ratio is 


given by 


R 

a 


= g<5 h T/L 

0(G) r 


Here, 6 T is a typical horizontal temperature difference over a typical 
length scale L and r is the lapse rate (T = -3T /3fc). We shall now try to 
obtain for baroclinic weather systems a lower bound of R . We know that V 
can hardly exceed the dry adiabatic lapse (g/C^) • Furthirmore, 

Max|G| = Hax(U/x, U^/L, fU) 


Here, U is a typical velocity associated with the scale L, f is the Coriolis 
parameter, and x is a typical time scale. Thus, overall 


(6 h T/L)C 

Min(R ) 

" max(U/x,U /L,fU) 


For large-scale flows in the middle latitudes, L ^ 106 m, U ^10 ms~l, 
f ^ 10~4 s-1, 0(G) ^ fU. Also, a modest estimate of the large scale 
temperature gradient in a baroclinic flow is 3 deg/1000 km; in addition. 

C = 1004 Jdeg kg , thus Min(R fl ) = 0(3). This means that in the 
l?h.s. of (11), the contribution of the terms associated with the horizontal 
temperature gradient typically dominate that associated with the vertical 
temperature gradient. The net result is 

g V H lnT * 3G/3z (11)* 

n v 

For the geostrophic case, G = (fv,-fu) and (11) 1 is recognized as an approxi- 
mate form of the thermal wind relation (e.g., HESS, 1959, p. 191). As long as 
significant baroclinicity exists, the approximation (11) ' continues to be valid 
for mesoscale flows with L ^ 105 m, U ^ 10 ms~l, | G | ^ U^/L and 6hT/L = 

0.3 deg/100 km. ^ 


Another useful form of (11) 


l.h.s. of (11) is actually gV^in^), 
horizontal gradient in x, y, p, t s' 


can be utilized if one recognizes that the 
Here the operator 


is so, note (HESS, 1959, pp. 260-264) that 


V„) is the 

ce (p is held constant)? To see why this 


V H T) Z = V H T) p + ( 3 T/3P)V h P 
Using hydrostatic (5) and the chain rule, we obtain 


3T/3P = -(1/Pg) 3T/3z 
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We also know from (10) that, (1/p) V P = G thus, overall, 

rl ^ 

SWz = sWp - W 3z 

The net result is 

g V H lnT) p = 3£/3z (11) " 

(11)" is an exact expression to calculate horizontal temperature gra- 
dients; nevertheless, to utilize it, its r.h.s. must be known at selected p 
levels. This requires knowledge of the pressure as a function of z. Typically, 
in the absence of rawinsonde, this is accomplished by utilizing a crude first 
guess of temperature from the radiometers, together with hydrostatic (5) and 
the equation of state (8). This results in a crude first guess of the pressure 
(typically + 10 mb). Equipped with this information, one can interpolate 
9G/3& which is observed in x, y, z, t space to an x, y, p, t coordinate. In 
practice, this is accomplished by interpolating to those z which correspond to 
constant p levels. Setting aside for the purpose of this discussion the 
standard errors associated with interpolations, there is an error associated 
with the fact that the pressure is inaccurately known; consequently, the z f s 
associated with the constant p levels are inaccurately known. 

We will now proceed to evaluate the above-mentioned errors. From 
hydrostatic we know that 

3 z/ Bln p = “RT v /g 

Integrating from sea level to some specified height (assuming for convenience 
p (sea level) ^ 1000 mb) we get 

z = (RT v /g)ln(p/1000) 

Here T is some vertically averaged temperature in the interval (o,z) and p 
is the v pressure at level z. In the lower troposphere ln(p/1000) = 0(1) and 
RT / g = 0(8 km) (see e.g., HESS, 1959, pp. 75-77). Assuming a worst case 
scenario that the errors in estimating T^ and T^ strictly from radiometric 
data are the same, we obtain for 6z (the error in z) that, 

6z = (R6T v /g)0(l) 

Taking 6T^ ^ 3°K we obtain 6z ^ 100 m. Now from (11)" and Taylor expansion 

6(3G/3z) % &z 3 2 G/ 3z 2 (12) 

'V V 

Furthermore, for Rossby number not too large from unity (essentially 
corresponding to large and mesoscale motion) g = 0(fU) thus, 

3 2 G/3z 2 = 0 (fU/D 2 ) (13) 

Here D is a typical vertical scale over which significant variation of 
9^u/ 3z^ are occurring. For the troposphere, D ^ 5 km, f 'v. 10“4 si, 

U v 10 ms“l, T Or 273 °K. Also as discussed above, Sz ^ 100 m. Substituting 
the above results in (12) and (13) and also taking into account (11)" we obtain 

6 V T) ^ 0(10~ 2 deg/100 km) (14) 

H p 

6V T) in (14) is an estimate of the error in the evaluation of the 
horizontal ^temperature gradient in p coordinates due to interpolation errors 
from x, y, z, t space to x, y, p, t space. As discussed above, these errors 
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are the result of our inaccurate knowledge of p. The error appears to be quite 
acceptable* Nevertheless, it should be remembered that our error estimate is 
quite sensitive to the choice of D, the vertical scale. At any rate, the 
algorithm to be described further utilizes (11) which is the exact form in x. 
y, z, t space rather than the approximate form (11)* or the form (11)" which is 
exact but requires interpolation to x, y, p, t space. 

Deduction of vertical virtual temperature gradients . So far we have 
shown how to find horizontal virtual temperature gradients. To find the 
vertical temperature gradients we will have to use the thermodynamic equation 
(6). A difficulty arises because while T v - T to within 1% (HESS, 1959, p. 

44) the contribution of the moisture to the horizontal virtual temperature 
gradient at the lower troposphere could be comparable to that of the horizontal 
temperature gradients. To overcome this difficulty we will now derive an 
alternative form of (6) containing only gradients of the virtual temperature. 

We start by noting that the continuity equation (1) and the equation of state 
(8) imply that 

- (1/ P)Dp/Dt = -RDT^/Dt + RT div* u 
This enables us to rewrite (6) as 

C DT/Dt - RDT /Dt + RT div*u = S, ( 15 ) 

p v v r\j n. 

From the definition of virtual temperature (9) we obtain 

DT/Dt = DT v /Dt - 0.61 qDT/Dt - 0.61 (Dq/Dt)T (16) 

New, under all meteorological conditions q ^ 10“2 and less. Thus, the second 
term of the r.h.s. of (16) is always negligible compared to the first term. 
Under conditions of strong moisture gradient (e.g., dry lines) the third term 
may be important and is therefore retained. However, in the third term, we may 
substitute T^ for T. Utilizing the above approximations we may substitute 
(16) in (15) V taking also into account the moisture equation (7) and the fact 

that C - C = R to obtain 
P v 

C V DT V / Dt - RT v divu - 0.61 C p S v T v = S fa (17) 

We next substitute (11) in (17) replacing the horizontal temperature gradient 
in (17) by (T v /g)9G/3z - (G/g)3T v /3z. The net result is 

C y (31nT v /3t + w31nT v /3z) = F - S h /T y < 18 ) 

Here F is a symbol for presumably observed quantities, i.e., 

F = l/g(u«3G/3z) + Rdivu - 0.61 CS, 

w is the modified vertical velocity given by 

w = w + (}fo-£)/g 

The horizontal velocity vector is denoted by [u^ = (u,v)] . 

The next step is to obtain explicit expressions for 31nT /3z which do 
not contain temperature tendencies. This is accomplished by applying the 
vector operator gV + £3/3z on both sides of (18). The result using (11) and 
calculus rules of the sort f'g = (gf ) 1 - g'f are terms like 

3/3t(gV H lnT v ) + £3/3t(3lnT v /9z) = (g v R + G <V 3z) 3lnT v / <*t 


(19) 

( 20 ) 
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gw 3/ 3z ( Vainly) + g (^jjW) SlnTy/ 3z = gV^(w31nT v ) / 3z 

Gw 3 / 3z ( 3lnT / 3z) + G( 3w/ 3z)31nT / 3z = G3/3z(w31nT / 3z) 

r\j * % V r U V 

Furthermore, from calculus 

G 3/ 3t ( 31nT /3z) = 3/3t(G31nT ,/3z) - (3G/3t) (3lnT /3z) 

Thus, 


(gV„ = G 3/ 3 z ) 31nT / 3t = 3/3t(gV + G3/3z)lnT 

H. Oj v n 'v v 

- (3£/3t) 31nT v /3z 

Using (11), 

(gV„ + G3/3z)31nT /3t = 3 2 G/(3z3t) - (3G/3t)31nT /3z 

H 'x v 'v ^ v 


Similarly, 


, „ 31nT „ 31nT 31nT n . 2 „ 31nT 

- « ( V>^ + gt r-rr- + *4 - 

3z 


Also utilizing (11) 

S. 1 

(g7 H 4G3/3z)— = — (gV H 4C3/3 Z )S h - (S^SG/Sz 
V V 

Overall, the net results are two separate estimates for the vertical 
temperature gradient, namely 


31nT 


C v (gV R w+£aw/3z - 9G/3t - w(3£/3z)-^ 
where H is given by 


- H 

a. 


( 21 ) 


H = (gV H +G3/3z)F - (S h /T v )3G/3z + (1/T v ) (gV H +G3/ 9z) S h 

Here F is given by (19). Since we have assumed that a first guess of T is 
available from the radiometers and is accurate to within + 3°K it is 
permissible to substitute this first guess in the r.h.s. of (22). The net 
result is that H is an observed vector function. 


( 22 ) 


It is now useful to put together the forms of the horizontal and vertical 
virtual temperature gradients and their dependence on the observed winds. They 
are : 


The generalized thermal wind relation (11) rewritten here as 

(g V u +G 3/ 3z) InT = 3G/9z 
n V 'V 

Equations for vertical temperature gradient (21) rewritten as 

A3lnT / 3z = H 

^ v ^ 

Here are the horizontal accelerations (with a minus sign), namely 

G = (-Du/Dt+fv+F 1 ; -Dv/Dt-fu+F 0 
'x 1 Z 

The term H is given by (22) and from (21) 

A = C v (gV H w+G3w/3z-3G/3t-w(3£/3z) 


(23) 

(24) 


(25) 


(26) 
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We close this section by noting that considerable simplifications of the 
expressions for A, H, and F may be realized if either the approximation (11) 1 
is used or pres sure' coordinates, are utilized. In that latter case (11)" may be 
employed. While we have shown that the above approximations are for the most 
part reasonable, we prefer the use of the more exact forms because the 
numerical solution of (23) and (24) does not become easier when (23) is 
replaced by (11) 1 or (11) 11 and approximate forms for A and H are utilized. 

Retrieval of the virtual temperature from the wind and radiances . 

Relations (23) and (24) contain information about spatial gradients of the 
virtual temperature. From this, as is discussed e.g., in GAL-CHEN (1978), a 
second-order three-dimensional Poisson-like partial differential equation for 
T may be obtained. To be solved in a limited domain, the boundary 
conditions (BC) need to be prescribed. GAL-CHEN (1978) has shown that 
Neumann- type boundary conditions (i.e., conditions on the virtual temperature 
gradient in the direction of the normal to the boundary) may be obtained from 
the observed wind (essentially from the components of £ and £ in the direction 
of the normal). Such a procedure appears to be better than using 
Dirichlet-type BC which require the specification of the virtual temperature 
itself on the boundaries. In the absence of radiosonde information, 
Dirichlet-type BC are usually known very crudely (either from radiometric data 
or from a guess from a larger scale model). 

Regardless of what type of BC are used, the use of (23) and (24) may not 
be optimal because it does not utilize the radiances from the infrared and 
microwave channels. Also, retrieval techniques based solely on (23) and (24) 
tacitly assume that the formulation of the dynamical equation (l)-(7) are 
infallible, i.e., that the retrieval errors of the virtual temperature would 
be attributed solely to observational uncertainties about the wind. 

We shall now proceed to develop a formulation which incorporates the 
observed radiances into the retrieval procedure. We start by noting that the 
radiative transfer equation may be reduced often to a Fredholm intergral 
equation of the first kind (e.g., WESTWATER and STRAND, 1972), namely, 

00 'Xj 

f B v (T)K(v, z )dz = I (27) 

Jo 

Here, v is the frequency, B v (T) is the flanck function, K are the weights and 
I v are observed radiances. The surface temperature contributions are included 
in the r.h.s. of (27). These contributions can be determined from the "window 
channel" measurements for space-borne radiometers and from the "big bang" 
cosmic background of 2.9°K for ground-based observations. For a well mixed 
gas, the function K(v, z ) is known except perhaps for a small temperature 
dependence. Traditional methods of determination of vertical temperature 
profile rely on solving (27) for various channels (frequencies) having 
different weights K(v, z ) . Thus, the contributions from different height layers 
can be varied and a degree of height resolution can be achieved. Extensive 
research (e.g., CHESTERS et al . , 1982) have demonstrated the limitations of 
such inversion techniques. In essence, the kernel K(v,z) acts as a vertical 
smoother (low pass filter). As a result, the retrieved temperature profile has 
a poor vertical resolution (at least in the troposphere). However, if (27) is 
combined with (23) and (24) the problem of vertical resolution is eliminated. 

In essence, the large vertical scales may be determined from (27) and the 
smaller vertical scales, which cannot be resolved by (27) would be determined 
from (23) and (24). 

Before we proceed with further mathematical developments of the idea 
outlined above, we note that (27) has been formulated for temperatures while 
(23) and (24) are valid for virtual temperature gradients. Furthermore, as has 
been noted before, the moisture contributions to the gradient may be 
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important. To express (27) in terms of the virtual temperature, observe that a 
Taylor expansion of B v (T v ) around T v taking into account (9) and the 
smallness of the virtual temperature correction would result in 

B v<V = B v (T) + (3B v / 3T)0.61 qT < 28 > 

Now, the second tern in the r.h. s. of (28) while small compared to the first, 
may not be neglected if we desire at least + 1°K accuracy for temperature 
retrievals; nevertheless, we may substitute the radiometers first guess about 
the moisture and temperature in the second term. To justify this approxi- 
mation, let us denote by <5 () the first guess retrieval errors. We may recall 
that 6q = + 6 gkg and <5T = + 3°K. We_also know that, T ^ 300°K and 
(for the lower troposphere) q ^ 10 gkg . Therefore, substituting the 
radiometers first guess in the second term of the r.h.s. of (28) would result 
in virtual temperature error ST^ of the order 

6T v = 0.61 q<5T + 061 T<$q 

Taking into account the order of magnitude of the various terms, the error is 
a£ most + 1°K. Furthermore, the contribution of this error to the radiances 
[I in (27)] is further reduced due to the averaging implied. Overall, we may 
substitute (28) in (27) approximating the second term in the r.h.s. of (28) by 
the radiometers first guess with the net result 
00 

[ B v (T v )K(v,z) = I v (29) 

. * 

where the moisture correction to the radiances I v have been absorbed in the 
term I v . 

The general retrieval algorithm may now be formulated as follows: Find 

a T v such that 

///[( g V H 4C9/3z)lnT v -3^/3z] 2 + (£31nT v /3z-g) 2 = Min (30a) 

subject to the constraint that 

00 

| B v (T v )K(v,z)dz = I v (30b) 

This is a familiar calculus of variation problem (C0URANT and HILBERT, 1953, 
Vol. 1, pp. 164-274) whose solution will not be discussed here. We note, 
however, that (30-a,b) attempts to satisfy the dynamical equations in the least 
square sense while enforcing the retrieved virtual temperature to satisfy 
everywhere the radiative transfer equation. 

A potential weakness of the retrieval algorithm is that the terms A and ^ 
involves calculating higher order derivative terms in both space and time. The 
estimate of such terms from the observed wind and its time history may be 
"noisy". To alleviate this problem one may use the Kalman filter approach 
(GHIL et al., 1980) where observations at more than two (or three) time levels 
are used to improve the estimate obtained from the solution of (30-a,b). 
Detailed examinations of the terms involved in (30-a,b) reveal that for the 
most part only two time levels are required. The calculation of A (equation 
26) requires knowledge of 9G/9t. Since ^ is acceleration, this requires 
knowledge of the wind at three time levels. Nevertheless, 9G/9t would be 
dropped out if we utilize pressure coordinates and relation (11)" or use the 
approximate form (ll) 1 . As noted earlier (11)" is exact but the use of pressure 
as a vertical coordinate requires some a priori knowledge of the pressure 
distribution. 
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